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ABSTRACT 

We present the results from several studies used to search for astrophysical sources of high-energy neutrinos 
using the Super- Kami okande-I (April 1996 to July 2001) neutrino-induced upward-going muon data. The data 
set consists of 2359 events with minimum energy 1.6 GeV, of which 1892 are through-going and 467 stop 
within the detector. The results of several independent analyses are presented, including searches for point 
sources using directional and temporal information and a search for signatures of cosmic-ray interactions with 
the interstellar medium in the upward-going muons. No statistically significant evidence for point sources or 
any diffuse flux from the plane of the galaxy was found, so specific limits on fluxes from likely point sources are 
calculated. The 90% C.L. upper limits on upward-going muon flux from astronomical sources which are located 
in the southern hemisphere and always under the horizon for Super-Kamiokande are (1 ~ 4) x lO^^^cm^^s^^. 
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1. INTRODUCTION 

There have been attempts to find extraterrestrial neutrino 
sources since the first underground experiments to observe 
cosmic-ray neutrinos weines et al. 1965, Achar et al.'i l965l) . 
Only two such sources have been observed so far. Me V elec- 
tron neutrinos from the Sun (Davis 2003; Koshiba 2003) and 
from SN1987A ( Hirata et al. 1987; Bionta et al 1987) . 

At GeV and higher energies, the observational method in- 
volves analyzing the arrival directions of neutrino-induced 
muons traversing underground detectors, since muons over 
about 1 GeV travel in roughly the same direction as their par- 
ent neutrinos. It is desirable to observe as high an energy 
sample of neutrinos as possible, as not only does this point- 
ing accuracy improve with energy, but astrophysical signals 
are expected to have harder spectra than atmospheric neutri- 
nos, gaining in signal-to-noise ratio in direct proportion to the 
energy threshold (Gaisser et al. 1995). Muons above a few 
GeV cannot be contained in the detector, so we restrict our 
search to only upward-going muons (UGMs) entering the de- 
tector. Such muons may be produced some distance away 
from the active detector volume, by charged current interac- 
tions of muon neutrinos (and anti-neutrinos) in the rock sur- 
rounding the detector. The effective volume for such instru- 
ments is thus roughly equal to the muon detection area times 
the range of the muon in the surrounding rock. The combina- 
tion of these factors means that any astrophysical point source 
of neutrinos found by this analysis would in the TeV energy 
range. 

Neutrinos may come from either astrophysical sources or 
cosmic ray interactions with Earth's atmosphere or the inter- 
stellar medium. Atmospheric neutrinos led to the observation 
of muon neutrino oscillations (Fukudaetal. 1998) but con- 
stitute an irreducible background in the search for astrophys- 
ical neutrinos. Furthermore, muons produced in the atmo- 
sphere penetrate from the surface at rates greater than muons 
from neutrino interactions in the surrounding rock, even in the 
deepest of mines. However, this background can be largely 
removed by considering only muons arriving from below the 
horizon, which can only be due to neutrino interactions. 

Given that cosmic rays are highly accelerated hadrons and 
neutrinos are produced as decay products of energetic hadron 
interactions, it follows that substantial neutrino fluxes must 
exist as well, although with orders of magnitude uncertainty in 
fluence. Many studies carried out over the last 25 years agree 
that a detector on the order of 10^ muon collection area 
would be required to see such sources ('Learned & MannheimI 
1^00) if they are at the lower end of the likely range of fluence. 
Thus, at the sensitivity level of Super-Kamiokande ("Super- 
K ") (about 1200 muon collecting area), the chance to un- 
ambiguously identify a source may be relatively small but cer- 
tainly not negligible. Additionally, bursting, highly beamed 
like objects such as gamma-ray bursts and blazars or other 
transient X-ray/7-ray sources for which the main mechanism 
is unknown and neutrino flux predictions are highly uncertain, 
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may produce fluxes of neutrinos detectable by Super-K . 

We report on analyses of the first five years of Super- 
K UGM data, searching for point sources and for diffuse flux 
from the galactic plane. No significant evidence for high en- 
ergy neutrino sources has been found. We report upper lim- 
its on neutrino-induced UGM fluxes from the directions of 
known high energy gamma ray sources, which in many cases 
are lower than those previously set. 

2. THE SUPER-KAMIOKANDE DETECTOR 

Super-K is a 50 kton ring imaging water Cherenkov detec- 
tor. Results presented here are derived from the "SK-I" run- 
ning period, which began when Super-K started taking data in 
April 1996 and ended when it was shut down for maintenance 
in July 2001. 

The detector is located in the Mozumi mine of the Kamioka 
Mining and Smelting Company, in Gifu prefecture, Japan, di- 
rectly underneath the peak of Mt. Ikenoyama with a minimum 
overburden of 2,700 m water equivalent. Its geodetic location 
is at 36.4° N, 137.3° E, and altitude 370 m. When looking 
at UGMs, Super-K 's northern location makes it sensitive to 
source locations on the celestial sphere with declination (5) 
less than 54°. 

The inner detector volume is a cylinder of radius 16.9 m and 
height 36.2 m, viewed during SK-I by 11,146 50 cm diame- 
ter inward-facing Hamamatsu photo-multiplier tubes (PMTs). 
Approximately 40% of the inner detector surface area is 
photo-cathode, with the remainder covered by light-absorbing 
black plastic. The outer detector volume, completely sur- 
rounding and optically isolated from the inner detector, is a 
layer of water 2.5 m thick viewed by 1,885 outward-facing 
20 cm PMTs. A more detailed description of the detector can 
be found in Fukuda et al. ( 2003). 

By using amplitude and timing information from the PMT 
signals, muon tracks in the inner detector can be reconstructed 
with tens of centimeters accuracy for endpoints and a few de- 
grees accuracy for direction, depending upon particle momen- 
tum (F uku da et al. 2003; Desai 2004; Ashie et al. 2005). The 
outer detector is used to identify particles entering and exiting 
the inner detector 

3. DATA AND BACKGROUND 
3.1. Data 

The dataset for these studies consists of all UGMs, both 
through-going and those that stopped within the inner detec- 
tor Contamination from cosmic ray muons prohibits the use 
of downward-going muons. Events were required to have 
path length in the detector greater than 7 m (corresponding to 
a muon energy greater than 1.6 GeV) to ensure good recon- 
struction accuracy. Parent neutrino energies, as determined by 
Monte Carlo simulation, are broadly distributed, with peaks 
around 10 GeV for stopping muons and 100 GeV for through- 
going muons. Further details on the data filtering and recon- 
struction may be found elsewhere (Ashie et al. 2005). 

Events with deposited energy greater than 1,750,000 equiv- 
alent photoelectrons (approximately 200 GeV deposited in 
the detector) were cut due to the likelihood of phototube 
saturati on. A separate stu dy has been done using those 
events dSwanson et alJ200^ . A dE/dx based method of tag- 
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ging higher energy UGMs is discussed in lDesail J2004h . but is 
not used in this paper's analysis. 

The total UGM livetime for the SK-I dataset was 1679.59 
days, about 87% of calendar time over the whole period. De- 
tector conditions were sufficiently stable throughout the live- 
time sample that we may assume them constant for this anal- 
ysis. Also, the exposure was sufficiently constant (x^ of 
75.3/95) in sidereal time to be taken as uniform for this pe- 
riod, as shown in Fig.Q 

The final SK-I UGM data set contains 1892 through-going 
muons and 467 stopping muons. 
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Fig. 1. — Sidereal time distribution of observed upward-going muons. 
There are no significant deviations from uniformity. 



3.2. Event Background 

In the energy range considered, the non-astrophysical back- 
ground to our studies is dominated by atmospheric neu- 
trinos. To understand if the Super-K data also contain 
some small contribution from astrophysical sources, source- 
free reference data are generated, used to test search meth- 
ods, and compared with the real data set. Two methods 
were used: Monte Carlo simulations of atmospheric neu- 
trino events, including detailed simulation of detector ef- 
fects; and the process of "bootstrapping" the real event 
data ( fSimpson & Maver-Hasse lwander 1986). 

In the bootstrap method, simulated event directions are gen- 
erated by randomly selecting uncorrected event times and di- 
rections from the real event database, and combining these 
time/direction pairs into a new simulated "event" which will 
project back to a different point on the sky than either of the 
real events from with the time and direction were drawn. This 
procedure is repeated until a statistically adequate sample is 
obtained. One thousand "fake" sky maps with the same to- 
tal number of UGMs as the database were generated to obtain 
the background estimates. Because real event times and direc- 
tions are used, any systematics due to UGM angular distribu- 
tion, detector angular acceptance, neutrino oscillation effects, 
live time variations between sky locations, etc., will be au- 
tomatically included in the background estimate, while any 
correlations present in the data will be destroyed. However, 
given the limited statistics of the real data in any given decli- 
nation band (shuffling in time only changes the event's right 



ascension), these backgrounds can be too granular, introduc- 
ing non-Poissonian effects. 

Analyses that did not employ the bootstrap method instead 
used Monte Carlo data to simulate the background. In or- 
der to compare expected background distributions with the 
2359 UGM events observed, 40 years of simulated atmo- 
spheric neutrino data (14610. Oequival ent live days) were gen- 
erated using Monte Carlo techniques jDesai et alJl2004^ . Be- 
cause the Monte Carlo events are generated without associ- 
ated times, random times within the same livetime range as 
the dataset were assigned to Monte Carlo events to generate 
fake data samples. 

4. TIME-INTEGRATED SEARCHES 

This section describes searches for possible astrophysical 
point sources of neutrino-induced UGMs looking just at the 
integrated data set, to see if there is a steady-state excess of 
events from some point in the sky. Previously the data were 
examined for evidence of such an excess from the centers of 
the Earth, S un, and Galaxy to set limits on possible WIMP 
annihilation (iD esai et alJ2004h . Time-dependent searches are 
described in §|5l 

4. 1 . Poisson Probability Sky Maps 

Several independent studies were performed on the SK-I 
dataset. The intuitive starting place for a point source search 
is a simple binned sky map. To generate such a sky map, 
the portion of the sky visible at Super-K for UGMs (i.e., be- 
low about 54° declination) is divided into bins, the number 
of events in each is then counted and a probability for finding 
the number of events due to background is generated for each 
bin. One study used 4° x 4° square bins which evenly tiled 
the v isible sky in an igloo pixelisation ( Crittenden & Turo^ 
11998ft . This size is somewhat sm aller than the more optimal 
4° half-angle cone discussed in § 14.31 but is the closest size 
igloo structure able to evenly map onto the sky. This square- 
bin study was performed four times, with bin center shifts of 
2° in right ascension and declination applied to account for 
potential sources near boundaries. 

Another study employed overlapping circular skymap bins, 
improving the sensitivity to potential sources near bin bound- 
aries and better matching the radial rather than rectangular 
expected shape of a source distribution. Though better able 
to avoid the potential effects of bin edges masking a real 
source, this oversampling by overlapping bins introduces non- 
Poissonian statistical effects which complicate error estima- 
tion. 

Both the square and circular bin studies used bootstrap es- 
timates of the background noise in each bin. They both pro- 
duced cumulative Poisson probability sky maps, an example 
of which is shown in Fig.|5]for one of the square bin searches. 
A point source would appear here as one or more bins with an 
anomalously high count relative to Poisson probability. How- 
ever, in both searches the data were consistent with the ex- 
pected background, as shown in the histogram of observed 
Poissonian probabilities in the lower half of the figure. A sam- 
ple drawn from a Poisson distribution should fit a power law 
of unit slope. The fit to the observed probability distribution 
has a slope of 1.06 ± 0.04, with a reduced — 0.83. No 
significantly improbable bins appeared in any of the sky maps 
described above. 

4.2. Clustering Analysis 
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Fig. 2. — Probability of the observed number of neutrinos in a skymap bin 
being a statistical fluctuation, for tiled 4° X 4° bins, shown on the sky (top) 
and as a histogram (bottom). The distribution of probabiHties (straight-line fit 
in bottom figure, slope = 1.06 it 0.04 and reduced = 0.83) is consistent 
with expectation for a Poisson distribution. 

Binned searches, such as the sky maps described above, 
treat the center of each bin as a test direction in which to 
look for a source, decreasing the sensitivity if a source re- 
ally lies near a bin boundary. Another approach is to look 
for statistically significant clustering of events on the sky map 
in a binning-free way. One such method (Ovama et al, 1989? 
lAmbrosio et al. 2001) involves using each event as a test di- 
rection. By examining the number of events in a cone around 
the test direction and comparing to the distribution produced 
by the background, one can look for regions of excess without 
boundary effects. 

This cluster-search algorithm was used on the Super- 
K UGM data. The expected background was calculated by 
dividing a Monte Carlo data sample equivalent to 40 years of 
Super-K operation into eight sets of 2359 simulated events to 
match the size of the actual data set. Using each event as a 
test direction, the number of events in a cone of half-angle 2° 
around each test direction were counted. The study was also 
done with an expanded cone of 3° . The distributions produced 
by the data show no significant deviation from those produced 
by the the average of the eight Monte Carlo-generated back- 
ground estimates. Fig. |3] The 2° search ha s a c omparable 
solid-angle to the binned all-sky surveys in §14.11 and the 3° 
matches the cones used for flux limits in § 14.31 

4.3. Specific Flux Limits 

None of the all-sky search methods described previously 
found a statistically significant excess of neutrino-induced 
muons over expectation from the atmospheric neutrino back- 
ground. We thus conclude that Super-K has not detected any 
sources of astrophysical neutrinos. 

Upper limits on the possible astrophysical neutrino-induced 
muon flux can be then be set for locations on the celestial 
sphere with 5 < 54°. Super-K 's northern location limits the 
exposure of the experiment to mainly sources in the Southern 
hemisphere. 

Flux limits have been calculated for a list of high-energy 




Number in Half-angle Cone 



Fig. 3. — Search for clustering around event directions. Plots show the 
number of occurrences of UGM clusters with multiplicity given on the hori- 
zontal axis. Upper plot is for 2° half-angle search cones, and lower plot is for 
3° half-angle cones. Data (points with error bars) are consistent with Monte 
Carlo expectation for background (histograms). 



astrophysical objects which have been described in the lit- 
erature as potential neutrino sources. For convenience of 
comparison, we start ed with the same hst of sources an- 
alyzed by MACRO (lAmbrosio et alJ 120011) . T o thes e we 
added the catalog of SGR's (discussed in § I5.1.2^ and 
other potential steady-state sources of TeV neutrinos such as 
LS5039 (Aharonianetal. 2005), WR20a (Bednarek 2005), 
IBS 1959H-650 (Reimeretal. 2005), B1509-58, B1706-44, 
and B1823-13 (>Link & Burgio^2005 ). For the flux estimates, 
UGMs within 3° of the source were selected as the signal. An 
angular window of 3° contains about 52% of the signal for 
an atmospheric neutrino spectrum and 92% of the signal for a 
neutrino spectrum. To evaluate the expected background, 
we used the atmospheric neutrino Monte-Carlo events. Each 
Monte Carlo event was assigned a time from that of the ob- 
served UGMs, in order to match the livetime distribution of 
the observed events. Then using this time and observed zenith 
and azimuthal angle, we can obtain the right ascension and 
declination for every Monte Carlo event, from which the an- 
gular separation between this event and any celestial object 
can be estimated. 

The numbers of signal and background were used to 
compute a 90% C.L. number of UGMs attributable to a 
background fluctuation using the procedure described in 
Caso et al. ( 1998), and that numbe r was used to compute a 
flux in cm^^s^^ as discussed in Fukuda et alJ in"999i) . The 
results are summarized in Tabled 

To minimize a priori issues, the approach of this paper has 
been to first do the all-sky search to see if anything anywhere 
is statistically significant (S I4.1> . Since there was not, UGM 
flux limits can be placed at specific locations of interest (Ta- 
ble[D. 

4.4. Diffuse Flux from Galactic Plane 

Galactic and extra-galactic cosmic rays can interact in with 
the interstellar medium (ISM) of our galaxy, producing a cas- 
cade of high energy particles, including neutrinos. Th ere have 
been various estimates of this flux in the literature JSteckeil 
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fT?79l iBe rezinskv et"an fT993t llngelman & ThunmarJ IT996t 
Eandia 20Q§. 

The density of hydrogen in our galaxy scales with height 
approximately as p oc with /iq = 0.26 kpc. Thus, 

most of the signal would come from the galactic plane. To 
search for UGMs from cosmic rays interacting with the galac- 
tic ISM, we first looked for an excess within ±10° around 
the Galactic plane. Since the flux of ISM-induced neutrinos 
is expected to be comparable to that of atmospheric neutri- 
nos only at energies above 100 GeV dlngelman & ThunmanI 
[T996ii) . we performed this search with only upward through- 
going muons. The galactic latitude distribution of upward 
through-going muons is shown in Fig.|4] Compared with the 
same distribution for the Monte Carlo data, there was no sta- 
tistically significant excess within ±10° of galactic latitude 
of zero degrees. Thus at the 90% CL, after subtraction of 
expected atmospheric neutrino background, less than 21 up- 
ward through-going muons observed at Super-K are from cos- 
mic rays interacting with the ISM. This flux limit is consis- 
tent with the predictions from the model by Berezinsk v et al.l 
lfT993 ). which when extrapolated to the size of Super-K pre- 
dicts about 2.3 events. We also note that recently a simi- 
lar search for high energy neutrinos from the galactic plane 
has b een performed with the AMANDA-II detector ( Kelley 
l20()5l) . 
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Fig. 4. — Galactic latitude distribution of upward through-going muon 
dataset compared with Monte Carlo expectation from atmospheric neutrinos. 
There is no significant excess near galactic latitude of 0° , which would be a 
signal of neutrinos being produced by cosmic ray interactions with the ISM. 



5. TIME-DEPENDENT SEARCHES 

In addition to direction, we also have timing information 
to aid in the search for point sources. This section describes 
searches for possible astrophysical point sources of neutrino- 
induced upward going muons which also consider the arrival 
times of the events in question, to be more sensitive to poten- 
tially transient sources than the "DC-style searches described 
in § 13 For these searches the timing cut used helps reduce 
the atmospheric neutrino background, allowing the use of a 
larger angular cone size (5° half-angle) as compared to time- 
integrated searches in order to maximize any possible signal. 



5.1. Triggered space-time correlation searches 

If a potential source is known to be transient or come in 
bursts from electromagnetic observations, the neutrino data 
can be examined in coincidence with these external "trigger" 
observations. The SK-I data set has aheady bee n searched for 
time coincidences with BATSE GRB triggers. jFukuda et alJ 

5.1.1. Markarian 501 Space-Time Correlation Study 

The flaring episode of Markarian 501 JWeekesll2000h be- 
tween February and October of 1997 was used to search 
for possible time coincidences with the SK-I UGM data. 
MRK 501 is a BL Lac object and a known TeV 7-ray source. 

A 5° half-angle cone centered on MRK 501 was used to 
count events recorded during periods of heightened activ- 
ity, which accounted for 13.5% of SK-I live time. Another 
67.5% was logged definitely before or after the flaring event. 
The remainder of the livetime corresponded to periods when 
MRK 501 was not observable by air Cherenkov telescopes, so 
its level of activity could not be determined and was not used 
in this analysis. 

One UGM was observed within this half-angle cone and 
during the flaring period. The estimated background which 
was estimated using the Monte-carlo data was 1.6. Given this 
background, at least seven events would be required in order 
to claim a 3cr signal detection. Therefore, the SK-I dataset 
provides no evidence for MRK 501 acting as a high-energy 
neutrino point source during its 1997 gamma ray outbursts. 

5.1.2. SGRs Space-Time Correlation Study 

Soft Gamma Repeaters (SGRs) are a small class (four con- 
firmed and three tentative candidates) of high energy tran- 
sient sources which emit bursts of X-rays and gamma-rays 
with durations lasting from 0. 1 to 500 sec, and are located i n 
the galactic plane (Hurley 2000; Woods & Thompson 2006^. 
Originally confused with GRBs, their burst properties differ 
from those of GRBs in two respects: SGR bursts are repet- 
itive, and have a softer spectral energy distribution. On rare 
occasions, SGRs have also emitted giant flares with hard spec- 
tra up to MeV energies. The total photon fluence received 
at earth varies from lO^^erg/cm^ to almost lO^-^erg/cm^. 
The number of bursts as a function of energy is given by 
^ « E^^-"^ . The currently favored model for burst emission 
from SGRs (Thompson & Duncan 1995) is magnetic field de- 
cay during crustquakes in neutron stars with B > lO^"* G 
(which are called "magnetar"). The flux of high-energy neu- 
trinos during a giant flare has aheady been estimated for the 
SGR 1806-20 flare on December 27, 2004, which was dur- 
ing SK-II (lokaetal. 2005; Halzenetal. 2005). There are 
also models for steady-state emission of high energy neutri- 
nos from selected SGRs ( Zhang et al.^2003) . 

We searched the SK-I data for space-time correlations be- 
tween UGMs and SGR bursts in an analysis similar to that 
used previously for GRBs (Fukudaet al. 2002). SGR burst 
times were obtained from the Interplanetary Network (IPN) 
catalog (K. Hurley, private communication). During SK- 
I only four SGRs were observed to emit bursts. These 
are SGR 1806-20 (33 bursts), SGR 1900H-14 (67 bursts), 
SGR 1627-41 (51 bursts), and SGR 1801-23 (2 bursts). We 
searched for coincidences within a window of AT = ±1 day 
and A9 = 5°. For this time window size, the total number of 
bursts from SGRs occurring in separate time bins was 74. The 
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space-time coincidences found between SGRs and UGMs are 
shown in Fig.|5] 
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Fig. 5. — Coincidences between SGRs and Super-K upward muons in 
a AT = ±1 day and AO = 30° window, where AT and A6» indicate 
the arrival time difference and angle respectively between the Super-K UGM 
and SGR. Each hatched box represents a correlation between an UGM and a 
group of SGR bursts which occur within a short time interval. The vertical 
mid-point of each box gives the direction between the observed UGM and 
the SGR, and the vertical size of the box is 2°, approximately the angular 
resolution for UGMs. Each symbol within a box gives the time of an SGR 
outburst relative to the UGM. The width of each box is the time interval 
between the first and last correlated SGR bursts. Within an angular search 
window of AS = 5° we found one UGM within hours of four different bursts 
from SGR 1900+14 (the bottom box) between 1998 May 26, 21:30:29 UT 
and 1998 May 27, 04:22:01 UT. 



No coincidences were found from SGR 1806-20, 
SGR 1627-41 and SGR 1801-23 and SGR 1627-41. 
One UGM was found within 2° of SGR 1900H-14 and 
within several hours of four bursts starting 1998 May 26, 
21:30:29 UT and ending 1998 May 27, 04:22:01 UT. The 
difference in times between the bursts and the UGM were 
-3.88, -3.75, 2.08, and 2.97 hours. 

The expected background in a AT = ±1 day window cal- 
culated using the Monte Carlo data is 0.013 events. Given 
that 74 different SGR bursts were examined, the total back- 
ground becomes 0.96. Therefore the UGM found coincident 
with SGR 1900-1-14 is consistent with the background and not 
statistically significant. 

5.2. Untriggered Searches for Burst Sources 

In the space-time correlation studies for UGMs described 
above, the reported times and sky coordinates of astrophysi- 
cal photon signals from other detectors were used as triggers 
to look for correlations in the SK-I neutrino data. We also 
looked for bursts of UGMs from astrophysical sources inde- 
pendent of photon observations, in order to detect transient 
sources of high energy neutrinos for which there may be no 
electromagnetic counterpart. The method used was to regard 
each observed UGM as a trigger itself, and look for other up- 
ward muon events arriving within a time window of one hour 
and an angular window of 5°. 

As we see from Fig. |6l using the entire SK-I dataset, one 
doublet was found to meet this criteria. Two UGMs separated 
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Fig. 6. — A scatter plot of the relative angle (A9) vs time-difference (A T) 
between consecutive UGMs (when arranged in chronological order) in SK-I, 
within a window of A T = 2 days and A6 = 20° . One UGM doublet was 
found in the signal window of A T = 1 hour and A9 = 5° . 



by 2.7° arrived within 19 minutes. The projected right as- 
cension and declination for these two events were 7''24'"0'', 
16°18' and 7''32'"24% 14°30' respectively. 

The arrival times of the events were Modified Julian Date 
= 50927.64412 and 50927.65757, or 1998 April 24, 15:27:32 
and 15:46:54 UT 

The expected rate of such coincidences is 2 x 10^^ per 
UGM. When multiplied by the number of trial factors (num- 
ber of UGMs minus one), the probability that this doublet 
comes from a background fluctuation is about 5%, which is 
not statistically significant. Thus we conclude that we do not 
see any evidence for bursts of UGMs from an astrophysical 
source in the SK-I dataset. 

6. SUMMARY 

Several independent methods for finding astrophysical 
sources of high energy neutrinos have been applied to the 
SK-I (1996 to 2001) upward-going muon data using both di- 
rectional and temporal information. None have uncovered a 
statistically significant excess above the atmospheric neutrino 
background as estimated from bootstrapped data and Monte 
Carlo simulations. Thus, upper limits were set on the upward- 
going muon flux. The 90% limits from astronomical sources 
which are located in southern hemisphere and always under 
the horizon for Super-K are (1 ~ 4) x 10" 
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TABLE 1 

The allsky surveys showed no significant point source, so these 
data are used to compute 90% flux upper limits within a 3° 
half-angle cone of the given equatorial coordinates of interest to 
astrophysical modeling. 



Source 


RA 


Dec 


# events 


Bkgnd 


90% c.l. p. flux limits 




(°) 


(°) 






(lO-i^cm-^s-i) 


SMC X-1 ' 


19.3 


-73.4 





3.0 


1.25 


T Mr X-2 1 


80.1 


-72.0 


5 


3.1 


3.49 


SN 1987A^ 


83.9 


-69.3 


2 


3.2 


1.9 


LMC X-4' 


83.2 


-66.4 


3 


3.3 


2.3 


GX301.2 ' 


186.7 


-62.8 


5 


3.5 


3.35 


Cen X-5 ' 


177.0 


-62.2 


4 


3.65 


2.73 


GX 304-1 ^ 


195.3 


-61.6 


4 


2.1 


3.33 


Cen X-3 ' 


170.3 


-60.5 





3.2 


1.27 


Cir X-1 ' 


230.2 


-57.2 


5 


3.3 


3.46 


21 J 1 fi17-^'? ' 


250.2 


-53.7 


3 


2.9 


1.99 


4U 1608-522 ' 


243.2 


-52.4 


5 


3.6 


3.10 


GX 339 4 ' 


255.7 


-48.8 


1 


2.4 


2.02 


Vela ' 


128.8 


-45.2 




2.7 


1.77 


GX 346-7 ' 


264.7 


-44.4 


4 


2.4 


4.29 


AR X-1 ' 


256.6 


-43.0 


2 


2.9 


2.66 




225 7 


-41 9 


3 


2 4 


3.65 


Vela X-1 ' 


135.5 


-40.6 


2 


1.7 


3.23 


2U 1700-37 ** 


226.0 


-37.8 





2.6 


1.90 


SGR X-4 ' 


275.9 


-30.4 


4 


2.0 


4.44 


LIO ^ 


269.4 


-30.1 


3 


2.0 


4.41 


GX l-H* ' 


263.0 


-24.7 


2 


1.5 


2.89 


SN 1604 ^ 


262.7 


-21.4 


3 


1.8 


3.79 


GX 9.9 ' 


262.9 


-17.0 


2 


1.6 


4.00 


Sco X-1 ' 


245.0 


-15.6 


3 


1.25 


5.43 


Aqr X-1 * 


310.0 


-0.9 





1.7 


2.5 


4U 336-1-01 1 


54.2 


0.6 


2 


1.5 


4.54 


Aql X-1 ' 


287.8 


0.6 





1.7 


2.53 


2U 1907-H02 


286.8 


2.3 


1 


1.9 


3.35 


Ser X-1 ' 


280.0 


5.0 


3 


1.6 


4.61 


SS433 ^ 


287.9 


5.0 


2 


1.5 


3.52 


2U 0613-H09 ' 

\J\J i. ^^yj y 


93.4 


9.1 





1.3 


2.7 




98.4 


17.8 





1.3 


2.94 


Crab ^ 


83.6 


22.0 





2.0 


3.09 


2U 035+30 ' 


58.8 


31.0 


3 


1.1 


6.78 




299.5 


35.2 


3 


1.6 


6.84 


Her X-1 ' 


254.4 


35.3 


2 


1.5 


6.94 


Mrk 421 * 


166.1 1 


38.2 


2 


1.0 


5.98 


Cvfi X-2 ' 


326.2 


38.3 


2 


0.8 


6.07 


Mrk 501 * 


253.5 


39.8 





1.0 


4.39 


Cyg X-3 ' 


308.1 


40.9 


1 


0.8 


6.66 


Per X-1 >° 


49.6 


41.5 





1.2 


4.66 


SGR 1806 " 


271.5 


-20.41 


1 


2.3 


2.82 


SGR 1900 ' ' 


286.8 


9.3 


1 


1.6 


3.62 


SGR 1627 " 


246.8 


-41.0 


2 


1.9 


3.11 


SGR 1801 " 


270.3 


-23.0 


1 


2.2 


2.78 


SGR 0525 " 


81.2 


-66.0 


5 


2.75 


3.64 


LS 5039 1 


276.5 


-14.8 


1 


1.2 


3.18 


WR20a' 


156 


-57.8 


3 


2.7 


2.51 


lES 1959-1-650 " 


300 


6.8 


3 


1.3 


6.31 


B 1509-58 3 


227.2 


-58 


6 


2.8 


3.67 


B 1706-44 3 


257.2 


-44.5 


1 


2.9 


2.14 


B1823-13 3 


275.6 


-13.2 


3 


1.35 


5.43 



'X-Ray Binary 
^SN Remnant 
'Pulsar 

''Emission Line Star 
'Dark Nebula 
^Cataclysmic binary 
^Symbiotic Star 
^Blazar 

'Wolf-Rayet Star 
'"Cluster 
"SGR 



